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^\ \ We present the theoretical phase diagrams of the classical Heisenberg fluid in an external mag- 

Q\ ■ netic field. A consistent account of correlations is carried out by the integral equation method. A 

nonmonotoneous effect of fields on the temperature of the gas-liquid critical point is found. Within 
the mean spherical approximation this nonmonotoneous behavior disappears for short-range enough 
spin-spin interactions. 
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The behavior of magnetic fluids in an external field compels more and more attention for the last years and has some 
peculiarities. In the presence of an external magnetic field the orientational (magnetic) phase transition is absent, 
but there are the first order transitions between ferromagnetic phases of different densities (e.g., gaseous and liquid 
phases). Physical properties of anisotropic fluids (to these belong also, besides magnetic fluids, nematic liquid crystals) 
are determined by the interplay between orientational and translational degrees of freedom. Therefore, varying the 
magnetic field it is possible to effect structural properties of magnetic fluids, in particular, to change the region of 
the gas-liquid coexistence. Such investigations with a calculation of phase diagrams were carried out for model spin 
systems within the mean field (MF) approximation It was found that for fluids of hard spheres carrying Ising 

spins an external magnetic field decreases the temperature of the gas-liquid critical point. On the other hand, the 
presence of isotropic van der Waals attractions between molecules can lead to the inverse effect . In Ref. the fluid 
of hard spheres with the classical Heisenberg spins and strong isotropic attractions was considered also. It was shown 
that at weak magnetic fields there can be two first order phase transitions in this model: gas-liquid and liquid-liquid. 
""O ■ In strong fields the weak liquid-liquid transition disappears. 

The need to account orientational-translational correlations for the description of physical properties of magnetic 
fluids stimulated studies of the external field effects by more complex techniques. The effect of an external field on 
the gas-liquid critical point was studied by the functional integration and Green function methods (|] for the quantum 
Heisenberg ferrofluid and by the Monte-Carlo and integral equation methods for the classical one The pair 

potentials of those models consisted of contributions of hard spheres and of the spin-spin interaction (the so-called 
ideal Heisenberg fluid). In these works the conclusion was made that an external magnetic field favors the phase 
■^j- ' separation, i.e., the application of the external field increases the gas- liquid critical temperature. Let us note that 
the results of Refs. ||,[| are obtained for quite strong fields. In our point of view, it was the effect of small fields 
that is of special interest. This follows from the fact that at small fields orientational fluctuations are large and the 
corresponding correlations have a long-range character. Therefore, small external influences can result in significant 
changes of macroscopic properties of magnetic fluids. Besides, the interplay between orientational and spatial ordering 
can lead to an interesting behavior of the gas-liquid critical point at small fields. In this letter we show that in the 
systems with the long-range enough spin-spin interactions the nonmonotoneous effect of the external fields on the 
gas-liquid critical temperature takes place. There is a temperature range, where weak external fields suppress the 
gas-liquid phase separation. By the integral equation method we show that while the range parameter of the model 
potential decreases this temperature interval first gets smaller and then disappears. 

We shall consider the model that was studied before for the case of zero external field 1^,0 . The pair potential 
(analogically to the papers H||) is a sum of the hard sphere potential ip(R\2) for spheres of diameter a and of the 
Heisenberg spin-spin potential $(i?i2, u)\, W2) 



^ ! *{R 12 ,W 1 ,U 2 ) = J(Rl2)&l-&2, (1) 

£ : J{R) = -K (za)2 *M-<R-<r)) : (2) 

za + 1 R/a 

where Si is a unit vector in the direction uii — (8i, ipi) of the magnetic dipole moment /x, referred to the uniform field 

Bo as the z direction. The potential of the particle interacting with the field is Vi = — fiBo cos Oi. In expression (|^) 
the coefficient <yZa ] , is chosen to make the integral 

ZCT+l ° 

I = ^ J duj-L J du 2 J dil 12 $(i? 12 ,Wi,W2)/K)/(a;2) (3) 

Rl2>T 
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(where f(u>) is a single-particle orientational distribution function) independent of za. The integral / describes a 
contribution of the spin-spin potential into the free energy functional within the MF approximation (see, for example, 
f|). Therefore, within the MF approach the model phase diagram is independent of za, if we use dimensionless 
units for the temperature t = k^T/K, the density rj — y^~, and the external field strength h — For the free 

energy of the hard sphere system we use the "quasi-exact" Carnahan-Starling expression || . The MF phase diagrams 
obtained by the well-known double-tangent construction are presented in Fig. |l|. One can see a nonmonotoneous effect 
of the external field on the temperature of the gas-liquid critical point at small values of h. With the increase of the 
external field strength the gas-liquid critical point first moves down (h — 0.1, 0.5), then moves up {h — 2, 5), and 
last (h — 10, 20, oo) grows higher than the top of the zero field binodal. It should be noted that the model potential 
in the limit za — ► belongs to a family of the so-called Kac potentials, for which the MF approach is accurate. 
Therefore, relying on the MF results (Fig. |l|) we can state that for long-range enough spin-spin interactions {za — > 0) 
the nonmonotoneous field effect on the critical temperature does take place. But for finite values of za we are forced 
to carry out more accurate investigation. 
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FIG. 1. Phase diagram of the Heisenberg fluid in the external magnetic field. MF approximation. The results are the same 
for any za. The thick lines constitute the phase diagram (the Curie line and the gas-liquid binodal) for h = 0. The thin lines 
are the gas-liquid coexistence lines for h — 0.1, 0.5, 2, 5, 10, 20, oo 



More consistent consideration of anisotropic fluids can be done on the basis of the integral equation method that 
allows to calculate both the one-particle distribution function and the pair distribution function. The task consists 
|iof in a solution of the anisotropic Ornstein-Zernike (OZ) equation 

h(l, 2) = c(l, 2) + J p(3)h(l, 3)c(3, 2)<f(3), (4) 

where d(3) = dR^duJs, p(l) = pf (u>i), h(l,2) and c(l,2) are the total and direct correlation functions of the system. 
Since Eq. (|J) contains the one-particle distribution function, we need (besides a closure relation for the anisotropic 
OZ equation) some additional relation for determination of p(l). It can be the first equation of the Bogolubov-Born- 
Green-Kirkwood-Yvon hierarchy (this was used in Refs. j|,[3| to obtain a numerical solution of Eq. ([I|)) or the Lovett 
equation for anisotropic fluids |llj] 

«(1) 



V W1 lnp(l) + V Wl ^U J c(l,2)V fc 



.P(2)d(2), 



(5) 



where V Wl is an angular gradient operator; v(l) is a potential of interaction with a uniform external field, its spherical 
harmonic expansion is of the form 



v(l) = -y]vtYio(tJi). 



(6) 



For a self-consistent solution of Eqs. (^) and ([}]) we shall use the method suggested in |Q,^2|. The method is based 
on the algebraic representation of the Lovett equation for uniaxial fluids. Using the general expansion for the direct 
correlation function of linear molecules 
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C(l, 2) = CnlW W^OO^A^*) (?) 

Til. 71 1 

,i„A 

and the exponential form of the one-particle distribution function 



/M = Z- 1 exp [5^y ro (a;)) 

\/>o / 



(8) 



we obtain following |7|jl3j an algebraic representation of the Lovett equation for a uniaxial fluid in the external field 

Ll = ^ ' Cl m Y mn L n + Vl — ^ ' ClraPm, + V h (9) 
mn m 

where all indices take the values greater or equal one, Vi = vi/k^T, C mn = J c]^ l0 {R)dR, Y mn = p (F m i(o;)l^ I 1 (a-')) w , 
(•••>« = / /(")(■■ -)^> ^1 = y/l{l + l)A h Vi = p v /Z(/ + l)(2/ + l)(P ; (cos0)) w , P>(cos0) is the Ith order Legendre 
polynomial. Let us note that the average values (Pi (cos 6))^ play the role of order parameters in anisotropic fluids. 
Relations (||) are accurate, and their use, as well as the use of the integro-differential equation (||) , does not introduce 
any approximation into the theory. 

Due to (0) it turns to be possible to obtain for the model (0)-(j^) an analytical solution of the anisotropic OZ 
equation (^) within the mean spherical closure 

c(l,2) = -$(i?i2,wi,w 2 )/fc B T 1 , i?i 2 >a, . . 

h(l,2) = -l, i?i2<CT. liUj 

Condition (jlTJ) for /i(l, 2) follows directly from the fact that hard spheres do not overlap. The mean spherical closure 
( |io| ) restricts correlation functions of our model to those of the form 

/(1,2)= £ fi lhm (Ri2)Y hm (u )l )Y l l m (u; 2 ), (11) 

(Zi 3 £2 = 0, 1), and representation ([)[) results in equalities 

A 1 = J 4 1 p(|y 11 (o;)| 2 ) w J Cni (R)dR+ 

A 1 = A 1 p{Y w (u)) u j c 111 (R)dR+-^ f . (12) 

Here we use the notations of Eq. ( pd|) for harmonics of the direct correlation function c(l,2); v% = pB^/^/Z. Thus, 
the use of the mean spherical closure yields for our model vanishing of coefficients A\ with I > 1 in Eq. (|8|) , and the 
self-consistent one-particle distribution function in the mean spherical approximation (MSA) takes the form 



f(uj) = exp (AxYioH) / / exp (A 1 Y 10 {u)) diu. (13) 

A uniaxial symmetry of our system leads to factorization of Eq. (|j) on the equations with different m. At m = ±1 

hu m (Rl2) = c llm (R 12 ) (14) 
+ P (\Yi m (u))\ 2 ^ J cii m (Ri 3 )hii m (R 32 )dR 3 . 

For m = we have a system of integral equations that after the Fourier transformation gains the matrix form 

Hij(k) = Cij(k) + J2 C »' Wpi'i'HMQi ( 15 ) 

i'j' 

where Hij(k) = hijo(k), Cij(k) = Cijo(k), ptj = p (Yio(u>)Yjo(u))) , indices take the values and 1. The problem 
of a self-consistent solution of the anisotropic OZ and Lovett equations has reduced at this stage to the solution of 
Eqs. (pj) and ( |l5| ) under conditions ( |l2"|) a nd the self-consistent f(w) given by relation (|l|). On the basis of the 
Wertheim-Baxter factorization method Wm one can find the analytical solution of such equations in the form of a 
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set of algebraic equations. The detailed derivation of similar solutions can be found in the literature, and therefore 
we omit any details and refer the reader to the previous publications 0,^2|- The explicit form of the solution is 
quite unwieldy and will be given in the extended version of this paper. Here we point out only that it is efficiently 
computable, and we use it for calculation of isotherms by the virial route to thermodynamics in order to locate the 
gas-liquid transition by the Maxwell construction. 




FIG. 2. Phase diagram of the Heisenberg fluid in the external field for za = 1 (a), za = 2 (b), za = 3 (c). Lines are results 
of the MSA. The thick line is the case of zero field. On the insets the vicinity of the ferrogas-ferroliquid critical point is shown, 
I and F mark isotropic and ferromagnetic domains. The thin lines are the gas-liquid binodals of the fluid in the magnetic field 
(the attached numbers are the values of the field h). Simulation data |j| (h = 0) for za = 1 (a) are shown as black circles (the 
gas-liquid coexistence) and diamonds (Curie points). It should be noted that ksT/K equals T*/6 from Ref. [[|. 

The configuration of the MSA zero-field phase diagrams slightly differs from that of the MF theory (Fig. [l]) . The 
MSA via the virial route to thermodynamics demonstrates (Fig. ||) the lack of the tricritical point JjJ in the Heisenberg 
fluid (|l])-(||) contrary to the MF and modified MF theories [g]. Within the MSA the Curie line joins the gas-liquid 
binodal at its vapor branch (see insets to Fig. |J). This result is in whole agreement with the available zero-field 
simulations for the same model |(|: the liquid phase is ferromagnetic and the gaseous phase is mainly paramagnetic, 
except in the neighborhood of the critical point, where the transition fcrroliquid-ferrogas takes place. The quantitative 
agreement with simulations is also quite perfect (see Fig. ||a). In the insets to the figure one can distinguish the critical 
point (CP, the top of the gas-liquid binodal) and the critical endpoint (CEP) in which the Curie line joins the binodal. 
In the temperature range from tcEP to tcp three spatially uniform phases (isotropic gas, ferrogas, ferroliquid) can 
exist. For long-range enough potentials this interval (tcp ~ ^cep) decreases with the decrease of za (see table and 
tends to zero for small values of za: the CP and the CEP coincide and form the tricritical point as is shown in Fig. [I]. 
In practice we can not distinguish the CP and the CEP by the MSA virial route already at za — 0.1. 

The effect of strong magnetic fields on the Heisenberg fluid (Q)-(|^) consists in a considerable increase of the critical 
temperature. Strong external fields spread the gas-liquid coexistence region on the phase diagram, in other words, 
favors the gas-liquid phase separation. This result agrees with the conclusions of Refs. But it follows from 

the MSA phase diagrams in Fig. || that small fields can suppress the gas-liquid transition at finite values of za (not 
only in the limit za — > 0). For example, one can see on Fig. ||a (za = 1) that the external field of strength h = 0.1 
totally removes the phase separation at temperatures from tcp(O.l) = 1.002 to icp(0) = 1.043. In the systems with 
more short-range anisotropic interactions this temperature interval decreases. For example, in Fig. |^b (za = 2) the 
temperature interval, where the external field h = 0.1 removes the gas-liquid separation, is much less. For short-range 
enough potentials even very weak fields do not suppress separation, e.g., for za = 3 the same field h = 0.1 increases 
tcp (Fig. ||c). Thus, the MSA predicts that the suppression effect of small fields on the gas- liquid separation gets weak 
and finally disappears for short-range enough spin-spin interactions. The verification of this conclusion we address to 
future investigations. 



za 






icp — icEP 


??CP — ??CEP 


2 


1.1913 


0.2119 


0.0016 


0.0185 


1 


1.0430 


0.1587 


0.0015 


0.0151 


0.5 


0.9811 


0.1316 


0.0008 


0.0063 


(MF theory) 


0.938 


0.117 









TABLE I. Coordinates of the gas-liquid critical point and its distance from the critical endpoint at different za 



4 



Acknowledgements. We are grateful to I. M. Mryglod for he compelled our attention to the problem. We thank 
J. J. Weis for sending us Ref. B prior to publication. We also thank H. Iro and R. Folk for giving us the paper Ref. j|] 
prior to publication and for useful discussion. 



E-mail: ccc@icmp.lviv.ua 
[1] T. Kawasaki, Progr. Theor. Phys. 58, 1357 (1977). 
[2] F. Schinagl, H. Iro, and R. Folk, to appear in Europ. Phys. Journ. B. 

[3] I. Vakarchuk, H. V. Ponedilok, and Yu. K. Rudavskii, Theor. Math. Phys. 58, 445 (1984). 
[4] F. Lado, E. Lomba, and J. J. Weis, to appear in Phys. Rev. E. 
[5] F. Lado and E. Lomba, Phys. Rev. Let. 80, 3535 (1998). 

[6] E. Lomba, J.-J. Weis, N. G. Almarza, F. Bresme, and G. Stell, Phys. Rev. E 49, 5169 (1994). 
[7] T. G. Sokolovska, Physica A 253, 459 (1998). 

[8] J. M. Tavares, M. M. Telo da Gama, P. I. C. Teixeira, J. J. Weis, and M. J. P. Nijmeijer, Phys. Rev. E 52, 1915 (1995). 

[9] N. F. Carnahan and K. E. Starling, J. Chem. Phys. 51, 635 (1969). 
[10] D. Henderson, in Fundamentals of Inhomogeneous Fluids (Dekker, New York, 1992), edited by D. Henderson, ch. 4. 
[11] R. A. Lovett, C. Y. Mou, and F. P. Buff, J. Chem. Phys. 65, 570 (1976). 

[12] M. F. Holovko and T. G. Sokolovska, Ukrainian Physical Journal 41, 933 (1996) (in Ukrainian). 
[13] T. G. Sokolovska and M. F. Holovko, Condensed Matter Physics 11 109 (1997). 
[14] M. S. Wertheim, J. Math. Phys. 5, 643 (1962). 



5 



